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Abstract 

The effect of fluctuation induced weakly first order phase transition known for 
three dimensional {3D) type I superconductors appears in a modified and strongly 
enhanced variant in thin (quasi-2D) superconducting films. The unusual thermody- 
namic properties of this new type of first order phase transitions and the possibility 
for an experimental verification of the effect are established and discussed. 



One of the outstanding problems of the phase transition theory of superconductivity is 
the so-called fluctuation induced weakly first order phase transition (WFOT) in three 
dimensional (3D) type I superconductors. It was shown |jl[ that the interaction between 
the fluctuations 61i of the magnetic field H and the order parameter ipl'x) may produce 
a fluctuation induced change of the order of the phase transition in 3D type I supercon- 
ductors: from a second order phase transition when the magnetic fluctuations 61i are 
neglected to a weak first order transition when the same magnetic fluctuations are taken 
into account and the mean magnetic field Hq = (H — 61i) is equal to zero. Note, that we 
consider nonmagnetic superconductors, where the magnetic induction B is equal to the 
magnetic field H. Besides, we discuss the phase transition point Tco corresponding to a 
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zero critical magnetic field, Hc{Tco) = 0, and our consideration will be not extended to 
the entire line [0 < Hc{T) < Hc{0)] of the superconducting phase transition. 

The same type of fluctuation induced first order phase transition was predicted also in 
other fields of physics, for example, in the theory of the early universe |0, ^ and on the basis 
of the de Gennes model of smectic A - nematic phase transitions in liquid crystals 0, ^ • It 
seems that under certain circumstances, i.e. suitable values of the characteristic lengths, 
such fluctuation induced first order transitions may occur in any system described by 
an Abelian-Higgs type of model 0, ^ H, where a scalar field like the superconducting 
order parameter ip^x) interacts with a massless vector gauge field. In a superconductor 
the massless gauge field is the vector potential A(x) of the magnetic field given for 3D 
systems by H(x) = x A and the Coulomb gauge divA = 0. A general description of the 
phase transition mechanism in Abelian-Higgs systems was given within the the scalar field 
electrodynamics The theoretical studies of this problem in superconductors have been 
performed with the help of the Ginzburg-Landau (GL) theory of superconductivity and the 
so-called dual model 0, |ll|), including the application of the renormalization group (for 
a review, see [§, |^)- The WFOT problem seems to be of general interest to a wide range 
of systems and, in particular, to recent aspects of superconductivity |jT2|, According 



to the theoretical predictions the effect in 3D superconductors is very weak |^ . For this 
reason an experimental search of the WFOT was attempted in smectic A liquid crystals 
by interface velocity measurements as well as by high-resolution heat capacity and X-ray 
experiments |T^. The results were interpreted [0, |15] in favour of a WFOT due to 



the interaction between the smectic scalar order parameter and the director vector of the 
nematic order, i.e. as a verification of the theoretically established fluctuation mechanism. 

We shall restrict our consideration to type I superconductors where the description of the 
WFOT can be given within the "approximation of spatially uniform order parameter" 
('0(x) ^ ip 1^). That is why we shall present a form of the GL free energy which 
can be used for the investigation of the WFOT in D-dimensional type I superconductors 
for a uniform order parameter ip. The 3D and 4D dimensionalities were analyzed in 
preceding works [jl], ^ and we shall focus on quasi-2D films. We shall demonstrate that 
the effect of WFOT is much stronger in films than in 3D samples and, therefore, could 
be experimentally observed (if it really exists). Besides, in quasi-2D superconductors the 
WFOT is described by a logarithmic term instead of a third order term in ip ^ , and 
this presents a special theoretical interest. 



In standard notations the GL free energy in case of a spatially uniform order parameter 
ip takes the form Fgl = (-^o + ^a) with 



Fo = V ^a\^P\' + "-\^\' j , (1) 

and 

PA = Y^J2jd''x [2p(V^)(5,, ^2 + {d,A, - d,A,)^] , (2) 
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where di = d/di, = 1,...,D; V = is the volume of the D- dimensional 

superconductor, H = |H|, a = ao{T — T^q) with ao(Tco) > 0, b(Tco) > 0, p = PolV'P 
with po = (Svre^/mc^). The quantities 2m and 2e are the effective mass and charge of 
the Cooper pairs respectively, and in superconductors of interest to the present inves- 
tigation 2m ^ 2me, where rrie is the electron rest mass. The concrete values of the 
initial (bare) critical temperature Tco and the Landau parameters ao and b given by the 
microscopic BCS theory |]16|, are not essential for our consideration within the general 
GL phenomenological approach. But we should keep in mind that the GL theory de- 
scribes only quasimacroscopic phenomena with characteristic lengths greater than the 
zero-temperature coherence length = (^^/4mQ;oTco)^^^; for a microscopic (BCS) justi- 
fication of this argument, see, e.g., . 



In the space of the wave vectors k, the upper cutoff A for = k is A ~ (l/^o) and 
can be extended to larger values only in case of special circumstances (calculation of 
cutoff independent integrals). Besides, the free energy (1) has additional restrictions 
which are also relevant to our consideration. Firstly, the approximation =const will 
be valid for well established type I superconductors, like Al, where the GL parameter k 
is much less than unity; k = A(T)/^(T). The coherence length ^(T) is given by ^(T) = 
(^o/I^K^ , where t = {T — Tco)/Tco, and the London penetration depth is A(T) = Ao/|t|^^^, 
where Aq = (mc^6/87re^aoT'co)^'^^ stands for the zero-temperature value of A. Secondly, 
we have to take into account the additional restriction on the GL theory for type I 
superconductors, namely, \T — Tco\ ^ k^Tco [jl6[, which comes from the general Landau 
condition \T — Tco\ <^ Tco and the requirement for a locality of the functional dependence 
between the supercurrent and the vector potential. The latter requirement is fulfilled 
for slow variations of the magnetic field and the vector potential at distances of order 
of ^0, i-e. for A(T) ^ ,^o- This means that the GL theory for well established type I 
superconductors with, for example, a GL parameter k ~ 10~^ -r- 10~^, is valid only for 
\t\ < 10~^ -7- 10~^K. The Ginzburg critical region for these superconductors is very small ( 
~ 10^^^ ^ lO^^^K [|, [1^) and can be safely ignored. All phenomena that can be reliably 
predicted within this theory are confined to a relatively narrow vicinity of the phase 
transition point (0 ~ 10"^^ < \t\ < 10"^ ^ iq""^). 

Thus we shall use the free energy (1) only for the description of phenomena in the tem- 
perature domain defined by the inequalities ^(T) ^ A(T) ^ which are fulfilled in the 
same vicinity of the phase transition point. Our heuristic arguments presented above have 
a reliable justification within the microscopic approach |]16 . 



Following preceding papers [||, ^ [T^ we shall integrate out the vector potential in order 
to obtain an effective free energy as a function of the uniform field ip. This integration 
is performed exactly within a simple one-loop expansion, as depicted in Fig.l, where the 
p— dependent term in eq. (2) is considered as a perturbation. The summation of the 
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infinite logarithmic series shown in Fig. 1 yields 

F^,^ = l{D-l)kBTj:^n[l + p{^)/e] , (3) 

k 

which substitutes Fa in Fgl- The total effective free energy density is given by /gg = 
{Fq + F^g)/V^. The function f^^i'ip) was investigated in details for D > 3 [|l], 0] when 
the effects of the term F^^ are small; a brief notice about the 2D case is also known. 

Let us focus the attention on the quasi-2D spatial dimensionality (3D thin films), where 
the film thickness Lq satisfies the condition A(T) ^ Lq ^ ^o- This condition is suitable 
for both theoretical and experimental investigations. Then our relatively thin supercon- 
ducting slab must obey the following conditions: 

e(T) » A(T) » Lo » ^0 . (4) 

It is easy to see that the k— summation in eq. (3) can be substituted with a 2D integration 
over the wave vector components parallel to the film surface and the transverse wave 
vector component is equal to zero. An upper cutoff A = (l/^o) is assumed for k = |k|. 
In contrast to the pure 2D case jl^ and studies in other spatial dimensionalities [|], |, |^ 



the result for the effective free energy of the quasi-2D superconductor is given by 

/eff = ^olV^P + f 1^1' - QI^Pln(Pol^lVA') , (5) 

where tq = (a + ksTpQ/ATcLQ), u = {b + kBTpl/A-nA'^LQ), and Q = {ksT pq/At:Lq). The 
evaluation of the contribution of the vector potential fluctuations, i.e. of the po-dependent 
terms in eq. (5) shows that these terms cannot be neglected for the entire temperature 
interval of validity of our consideration. Besides, the evaluation of the new terms in the 
free energy is useful as a demonstration of the important role of the cutoff A and the 
necessity of the consistent choice A ^ (l/^o)- Another choice of A immediately leads to 
wrong predictions. 

The free energy /(y?) = [/gff(v5)/'u] with ip = > Q is depicted in Fig. 2 for several 
values of r = [tq + Q\yi{K^ / po) /u], and q = {Q /u). Figure 2 shows a well estabhshed first 
order phase transition. We shall briefly summarize the analytical treatment of the free 
energy (5) in terms of the simple notations /, r, and q. The analysis is similar to that 
for first order phase transitions described by a third order )term; for details, see p. 
Let us denote the solutions of the equation of state {df /dip) = by v^o(= |^o|) - the 
possible phases. The normal phase solution (y^o = 0) exists and gives a minimum of / for 
all T > 0. Moreover, the necessary stability condition f"{(po) = {d"^ f /d'^'^)-fio > exhibits 
a singular behaviour /"(v^o 0) — > {+oo), i.e. a logarithmic divergence. This peculiar 
property does not create troubles about the physical meaning of the effective free energy. 
Rather this result shows that in thin films the normal state can occur at any T > as a 
stable phase above the equilibrium phase transition temperature Tc ^ Tco, or, if certain 
experimental conditions are satisfied, as a metastable phase for < T < T^. 
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The superconducting (Meissner) phase y^o > is described by the equation 

{r-q) + ^l- q\n^l = . (6) 

The solution of this equation, i.e. the equihbrium order parameter (/jq > is shown in 
Fig. 3 for values of the parameters r and q which allow the existence of a first order 
transition. The check of the necessary {f'\ipo) > 0) and the sufficient (/(v^o) < 0) 
stability conditions for y^o > show that the former is fulfilled for ^Jq > g, whereas the 
latter condition is fulfilled for (p^ > 2q. Therefore, the Meissner phase can be overheated 
above the equilibrium transition temperature Tc of the phase transition defined by the 
equation v^o(^c) = '^<l{Tc)- The equihbrium order parameter jump is equal to (2gc)^^^, 
where qc = q{Tc)- This overheating may continue up to the temperature defined by 
(pl{Tci) = q{Tci) , where the order parameter becomes equal to - the lowest nonzero 
value of the order parameter, for which it gives a minimum of /. These features of the 
order parameter are shown in Fig. 3, where the "metastability extension" of the order 
parameter profile is also drawn for a suitable choice of the parameter q. Note, that one 
can easily show in an analytical way that stable solutions < y^o < 1 of the eq. (6) always 
exist, at least, for r < 0. 

The latent heat L{Tc) = TcAS{Tc) and the specific heat jump AC(Tc) = Tc{dS/dT)T^ at 
the equilibrium phase transition point Tc can be calculated from the free energy (5). We 
shall be interested in the temperature size of this first order transition, namely, in the 
ratio (AT)c = \L(Tc)/C(Tc)\. Taking into account only the temperature dependence of 
the effective free energy /gg on the Landau parameter a, the temperature size (AT)c can 
be correctly evaluated. Our result is 

(AT), = ^ , (7) 
and after the substitution ao = {h^ / AniTcQS,'^)^^'^ and Q = [ksT p^/ At:Is?Lq) we have 

When we put in the above expression the numerical values of the fundamental constants 
and the tabulated data for Al (T, = 1.19 K, and ^(0) = 1.6 x 10""^ cm [|]), we obtain 
that (AT)c ~ 10~^/To. Note, that our approach is valid, if the thickness Lq is greater 
than ^0 and less than A(T). A choice of Lq consistent with all theoretical requirements 
is To ~ 10 X ^0, or, for Al, Lq ~ 10"^ cm, which yields (AT)c ~ 10"^ K. As the usual 
GL theory for Al is restricted within a temperature interval of size 10~^ K, our result 
indicates that the size of the WFOT covers totally the domain of validity of our unusual 
theory, namely we are faced with a normal size FOT rather than with a WFOT. Therefore, 
it should be emphasized that if the gauge mechanism of a fluctuation induction of the 
phase transitions order change really exists, this should be seen in suitable experiments 
with pure samples of good type I superconducting fllms of thickness about lO/zm. The 
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same result for the temperature size of the FOT can be obtained by the estimation of the 



parameter ao given by the BCS theory |16|. 
In conclusion, we have three notes: 

(1) Our approach to quasi-2D superconductors allows to use the 3D values of the original 
GL parameters ao and b and keeps the consideration far from dangerous effects of distruc- 
tion of the order by 2D fluctuations. In this respect as well as in the values of the effective 
parameters vq, u, and Q, our effective free energy (5) is quite different from the pure 2D 
free energy known from a preceding paper of Lovesey fl^. That is why we are able to 



make the prediction about the enhancement of the WFOT in quasi-2D superconductors 
(or, which is the same, "almost 3D" superconductors). We should stress that for pure 
2D case the Landau parameters and b are different from those corresponding to 3D 
superconductors and, hence, the evaluation of the same FOT effect in two dimensions 
needs a consistent treatment of these parameters. Moreover, the investigation of very 
thin (almost 2D) films requires an evaluation of the temperature region where the order 
parameter fluctuations will destroy the phase transition. 

For these reasons our results cannot be extended in a straightforward way to very thin 
films by the respective decrease of the film thickness Lq. Within the phenomenological 



approach to dimensional crossover phenomena , the 3D Landau parameters are related 
to the 2D ones by X(2D) = liKiLo->o(-^-^o)) where X = {ao,b). This consideration should 
be compared with the BCS predictions for the 2D values of the free energy parameters. 

(2) We have followed a method of integration of the magnetic fluctuations proposed in 
preceding papers ||T], , in which the order parameter fluctuations are not totally ignored 
and the "photon mass" p{ip) generated by the effect of the Higgs field ip on the vector 
potential A is not an equilibrium quantity. The special feature of the present method 
is that the uniform fluctuation Sip of the uniform "field" ip interacts with the magnetic 
fluctuations and this interaction does ensure the fluctuation mechanism of the FOT tran- 
sition discussed in the present and all preceding works. The WFOT in higher spatial 
dimensionalities and the FOT in our case are caused by a high-order fluctuation interac- 
tion [{6il!)'^{6Aj)'^] between the order parameter and the magnetic field fluctuations rather 
than the magnetic fluctuations alone. 

(3) The theoretical arguments in favour of the FOT within the approximation ip = const 
below the equilibrium transition point are not precisely the same as those obtained by 
various renormalization group investigations in one loop |^, |l^ and higher-order [^, 

approximations in the loop expansion, where another type of singularities of the 
perturbation series for the GL free energy are relevant. Within the present approximation 
for ijj the divergent terms, which are usually relevant for the renormalization group studies 
and are present in all diagrammes in Fig. 1, have been summed up to a finite sum. 
Therefore, the arguments in favour (or against) a WFOT within the approximation of 
uniform scalar field ip cannot be used to justify or reject renormalization group predictions 
above the equilibrium transition point. 
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Figure captions 



Fig. 1. One loop logarithmic series of diagrammes: the black squares denote the vertex 
part of the A^— term, and the sohd lines denote the correlation function < |^j(k)p >. 

Fig. 2. The shape of the effective free energy f ior q — 10~^: hne (1) corresponds to 
r = -9.3 X 10-^ (2) - to r = -9.4 x 10-^ (3) - to r = -9.515 x 10-^ and (4) - to 
r = -9.7 X lO"! 

Fig. 3. The order parameter profile for q — 10^^. The black squares indicate the stable 
superconducting states and the white squares indicate metastable superconducting states 
above the equilibrium transition point. 
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